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Abstract. We study the coherent dynamics of an excited two-level atom in a 
vacuum optomechanical cavity and find that the original atom-cavity Rabi oscillation 
is sinusoidally modulated by the light-mechanics coupling as the Rabi splitting is on 
resonance with the mechanical mode. We develop an analytic model in a three- 
dimensional Hilbert subspace to explain this phenomenon and employ numerical 
simulations of the density-matrix master equation to confirm our analysis. We also 
show that the modulated Rabi oscillation survives in presence of dissipations and other 
non-ideal factors. 


1. Introduction 

Optomechanics and cavity quantum electrodynamics (QED) both study light-matter 
interactions. While cavity QED involves light interacting with nearly-resonant systems 
such as atoms Dia. optomechanics treats highly off-resonant couplings between light 
and mechanical objects [SI HJ [5l El [7] . The rapid developments in optomechanics towards 
the strong-coupling regime make it possible to explore the quantum effects of single¬ 
photons on macroscopic resonators [SI El cni [11112111 US], and hybrid systems consisting 
of optomechanics and cavity QED components are expected to exhibit rich and novel 
features in the quantum regime [H [161 [13 [IH]- 

As is well known, an excited two-level atom in a vacuum cavity undergoes Rabi 
oscillation at a frequency proportional to the atom-photon coupling. This is one excellent 
example illustrating how the cavity modifies the optical properties of the atom, since in 
free space the atom would simply decay exponentially. Adding a mechanical resonator 
to the atom-cavity system via optomechanical coupling further changes the atomic 
behavior, even though the mechanical resonator does not directly interact with the atom. 
In this paper, we study the coherent dynamics of an excited two-level atom in a vacuum 
optomechanical cavity. We find that as the atom-cavity Rabi splitting is on resonance 
with the mechanical mode, the original Rabi oscillation is sinusoidally modulated by the 
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optomechanical (light-mechanics) interaction. This sinusoidal modulation of the Rabi 
oscillation is non-trivial since it occurs in absence of any initial photons and phonons. An 
analytic model in a three-dimensional Hilbert subspace of the hybrid system is developed 
to explain this intriguing phenomenon, and numerical simulations of the density-matrix 
master equation are employed to conhrm the analytical analysis and to demonstrate 
that the result qualitatively holds even as dissipations and other non-ideal factors are 
taken into account. 


mechanical 

resonator 



Figure 1. Schematic of a two-level atom inside an optomechanical cavity. 


2. Model and Hamiltonian 


We consider the setup in Fig. 1 , where a two-level atom is placed in a cavity with 
a movable end mirror as the mechanical resonator. The single optical cavity mode is 
coupled to both the atom (via dipole-moment interaction) and the mechanical resonator 
(via radiation pressure). The hybrid atom-optomechanical system is described by the 
Hamiltonian (h = 1) |T6] 


H = Ucc’c + -h Qca (d+c + d_c'^) + Umb^b - gcmc’c (b -f 6'^) 


( 1 ) 


Here c and b are the annihilation operators for the optical held (of frequency ujc) and 
the mechanical resonator (of frequency cOm), respectively, d^ = |e)„ „ (e| — \g)aa{ 9 \y 
d+ = dl = |e)^ ^ ((7I, with |e)^ (IS')a) ^^e excited (ground) state of the atom, Ua is the 
transition frequency between the atomic states, and gem (fi'ca) characterizes the strength 
of the light-mechanics (light-atom) interaction. In this letter, we use subscripts “c”, 
“a”, “m” to respectively denote “optical (cavity)”, “atomic” and “mechanical” states 
or parameters. The dimension of the hybrid-system Hilbert space is inhnite. However, 
as we will show later, if the system is initially in |0)^ |e)^ |0)m! i-®-; with the atom fully 
excited while the optical and the mechnical modes in the ground states, then during the 
evolution the system is conhned to a three-dimensional Hilbert subspace, which allows 
for an analytical description of the coherent dynamics. 

We note that the photon number plus the population in the atomic excited state is 
a constant of motion. Thus, starting from | 0 )^ |e)^ | 0 )^ , the quantum state of the hybrid 
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system at any time t stays in the Hilbert subspace \g)^ , |0)^ where |1)^ 

is the one-photon Fock state, “{}” stands for the subspace spanned by the basis vectors 
inside, and 'Hm represents the Hilbert space of the mechanical mode. Assuming that 
the atom and the cavity are on resonance, i.e., Ua = ujc, we construct the eigenstates of 
the atom-cavity part of the Hamiltonian [Hca = ooc&c + u:aC^z/‘^ + Qca (d+c - 1 - d_cf)] in 
the two-dimensional subspace {|1)^ \g)^ , |0)^ |e)^} as 

I , ^ _ \^)c\9)a± |0 )c|e)a 

\±)ca ^ 

which correspond to the eigenvalues of ujal‘^^gca- In this new basis Hca is diagonalized, 

i.e., 

Hca = Y + (3) 

with d; = 1+),, ,,(+! - |-)eaca(-|- Moreover, since ctc|±)^„ = (|+),, + h) J/2, 
one has cic = (^1 -|- d^ -|- d^) /2, where d^ = d^ = |+)ca ca (“I- Thus in the subspace 
|0)^ |e)^} 0 Hmi the total Hamiltonian ([T]) becomes H = a;a/2 -|- f^cad' -f 
ujrnb^b — gem (l + d^ -|- d^^ (b + /2. By introducing the displaced phonon operator 

b' = b — gem! the Hamiltonian is re-written as 

H = gcaS'^ + oOmb'^^y - ^ (d+ a'_\ (h' + b'A - (d+ a'_) , 

Z ZijJyyi 

where a constant energy term ujal‘2. — glm! has been discarded. Further assuming 

that gem ^ ^m and Um ~ ‘^gca) we neglect the term proportional to g^.^ and make the 
rotating wave approximation in the above Hamiltonian to obtain H Heff, with 

Heff = gcaS'^ + UJmb'^y “ ^ (d+&' , ^4) 

This Hamiltonian describes a two-level “polariton” coupled to a bosonic mode, analogous 
to the original atom-cavity Hamiltonian Hca- 

3. Sinusoidally Modulated Rabi Oscillation 

In the previous section, we derived an effective Hamiltonian Heff in dl]) for the 
hybrid system where initially the atom is fully excited and both the optical and the 
mechanical modes are in their ground states. At hrst glance, since Re// is a J-C 
type Hamiltonian, it might suggest a simple Rabi oscillation of frequency gem between 
two basis vectors, analogous to that in a regular atom-cavity system (i.e., without 
the mechanical resonator). This would be true if the hybrid system were initially in 
l+)ca |0)mo with |0)^, being the ground state of the displaced mechanical mode b'. 
However, the initial state of the system is not |+)ca |0)m'! rather lO)^ |e)^ |0)^. 
As will become clear below, starting from |0)^ |e)^ |0)m> evolution of the hybrid 

system involves three (rather than two) basis vectors, which in turn leads to a non¬ 
trivial modulation of the original atom-cavity Rabi oscillation by the optomechanical 
coupling. 
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We first take a deeper look at the initial state |0)^ |e)^ |0)m- Although rigorously the 
ground state |0)^ of the original mechanical mode b is the coherent state of the displaced 
mode b' with the eigenvalue —gem/ we can well approximate |0)^ ~ |0)^, since it 

has been assumed that gem ^ ^m- According to ([2]), |0)^ |e)^ is a superposition of |+)ca 
and |—)ca, and thus the initial state |0)^ |e)^ |0)^ ~ |0)^ |e)^ |0)^, is a superposition of 
l+)ca |0)m' \~)ca |0)m'- Further noting that |-)^^ |0)^, is an eigenstate of -ffe//, and 

Heff couples |+)ca |0)^/ only to |—)ca |l)m' ^Fe versa, we conclude that the coherent 
dynamics of the system is constrained to the three-dimensional Hilbert subspace of 

= {|-)„ |0>.„.. l+>„ |0>.„.. |->o. |1>„J . (5) 

which allows us to obtain an analytic solution for the evolution of the hybrid system. 

To this end, we diagonalize Heff in the subspace 


Heff = 


~gca 


0)(0 


-l- f 2 _l 


+ + 


-h Vt_ 


( 6 ) 


where 


0 ) =|-)ca| 0 ). 


+) - AA(|-)c„ |l)m/ - g | 0 )^,), 

-) =Af{g |-)^^ |1)^, -f |0)^,), 

and i ^9) ^9 \/‘^gi^ T gcmi V [^9 (n^m ‘^gca)\ /gc 


(7) 
= 

(1 -|- g‘^) With the initial state of the system decomposed [by inverting ([2]) and ([7])] 
in terms of the three eigenvectors of Heff as 


M (g 

+)- 

-)) + 

|e)JOL, =-^ 




we get the quantum state at any time t 

I / A^ . i^m + ng , J\f _ 

|V.(t)) = -^r,e . |+) + ^e 


V2 




To calculate the probability for the atom to stay in the excited state, we convert 


±y and Oy in 1-0 (f)) back to the “original” basis vectors of |0)^ |e)^ |0)m'’ |7)c \9)a |0)mo 
0)c \^)a |7)m' |7)c \9)a \'^)m'i which cau be done by simply substituting ([2]) into ([7]). 

Here we do not present the detailed expression of \'4) {t)) in the original basis, but only 
write down its probability amplitudes Aoei for |0)^ |e)^ |1)^, and ^oeo for |0)^ |e)^ |0)m'- 






oei = - iM^ge 


e 2 


-AoeO — 


+ e' 


'i'Qoa^ 


The probability Pg to hnd the atom in the excited state is given by |^oei|^ + I'Aoeol^- In 
terms of the more widely used poulation inversion AP = Pg — Pg, one has 

AP = cos iyteat) cos ^ ^ siu (Hcaf) siu 


( 8 ) 
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where 


^ca 


UJr. 


+ 9cai i — 


1 — 7f 

l + yf 


and Pg = 1 — Pe V the probability for the atom to be in the gronnd state. Eq. ([8]) is 
drastically different from the standard Rabi oscillation APca = cos {2gcat) in the atom- 
cavity system governed by Hamiltonian Hca [©I- Particnlarly, if the atom-cavity Rabi 
splitting is on resonance with the mechanical mode, i.e., Um = ^.gca, we can simplify (IHll 
into 


AF = cos {2gcat) cos , (9) 

which shows that the optomechanical coupling modulates the original atom-cavity Rabi 
oscilation with a cosine envelope function. This sinusoidally modulated Rabi oscillation 
is our central result, and its derivation has involved a fairly large amount of mathematics, 
so in next paragraph we will provide a more physically intuitive argument. 

To do this, we analyze the action of Fe// in ® on the three dimensional 
Hilbert subspace in ©. We note that Heff induces two distinct oscillatory 

transitions: on a short time scale of gff^, the atom-cavity interaction {gcaSf) yields the 
Rabi oscillation between |0)Je)„ |0)^, = (|+)^„ - \-)ca) \^)m> /and |l)jc/)j0)^, = 
(l+)ca + l“)ca) |0)m' /or equivalently, the oscillation of the dynamical phase 
difference between the |+)ca |0)m' and |—)ca|0)m' components; on a long time scale 
of the optomechanical coupling [—gem /2] causes the “polariton- 

phonon” oscillation between |+)ca |9)m' and |—At t = 0, the system is in 
|0)c \^)a |9)m' = (l+)ca “ l“)ca) \^) m' / and the atom-cavity interaction generates the 
Rabi oscillation on the short time scale. After a duration of 'k jgem^ l+)ca |9)m' evolves 
into |—)ca |l)m' "^ae to the optomechanical coupling, and the state of the system becomes 

is a phase difference between the two Fock states 
of the (displaced) mechanical mode. Around this time {t ~ T^/gcm) there is no atom- 
cavity Rabi oscillation since the |+)(,q |0)^, component dose not exist in the system 
state. After another duration of Tr/gcm, he., at f = 2n/gem, l“)ca l^)m' evolves back to 
l+)ca |9)m'! atom-cavity Rabi oscillation re-appears. For other t G {0,271 /gem), 

partial Rabi oscillations (i.e., without completely reaching the ground or the excited 
states of the atom) occur. This pattern repeats in every period of 2n/gem, leading 
to the periodically-modulated Rabi oscillation in ©. We notice that the envelope 
function in © has a period of Att/ gem, but the Rabi oscillation on the short time scale 
does not distinguish positive and negative values of the envelope function, and thus the 
“modulation period” is indeed 277/gem, agreeing with our analysis above. We further 
remark that the involvement of three (rather than two) basis vectors in is crucial to 
the sinusoidally modulated Rabi oscillation, because it is the optomechanical transition 
between |+)ca |0)m' l“)cal^)m' gives rise to the modulation of the regular 

atom-cavity Rabi oscillation in the subspace of {|+)ca |0)m' > l“)ca 
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Figure 2. Temporal evolution of (a) the atomic population inversion, and (b) the 
mean phonon number, with gem = OTw^, Wa = Wc, gca = O.Swm, and no thermal 
effects. The three curves in (a) are almost indistinguishable from each other, and the 
sinusoidally modulated Rabi oscillation is nearly perfect. 


4. Numerical Simulations and Non-Ideal Situations 


Eq. (|9]), which illustrates the sinusoidal modulation of the Rabi oscillation in the 
atom-optomechanical system, has been derived with the following assumptions: no 
dissipations, Uc = coa, oJm = ‘^dca, Qcm -C oJmi and deviations from these ideal conditions 
are expected to affect the dynamics. For example, if Um is slightly detuned from 2gca 
(and thus ^ 7 ^ 0), then AP is determined by Eq. (|H]) [rather than Eq. (|9])], where the 
second term on the right-hand side perturbs the “ideal” dynamics manifested in the 
first term. To more comprehensively understand the robustness of the modulated Rabi 
oscillation, we resort to numerical simulations of the density-matrix (p) master equation 
[19] for the hybrid system 
dp 1 


dt 


i 


H,p -KL[d]p-'jL[a-]p-{nth + l) pL b p-uthpL 6 ^ p, (10) 


where H is the unapproximated Hamiltonian in ([T|) , L [d] p = d^dp/2 — dpo^ + pd^d/2, k, 
7 , p are the decay rates of the optical cavity, the atom, and the mechanical resonator, 
respectively, and nth is fhe thermal occupation of the mechanical heat-bath mode at 
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Figure 3. Temporal evolution of (a) the atomic population inversion, and (b) the 
mean phonon number, with gem = O.lwm, Wa = Wc, gca = O.dSwm, and no thermal 
effects. In (a), the curve for Eq. (jH]) essentially coincides with the numerical one, which 
shows a non-ideal modulated Rabi oscillation. 

frequency LUm- It has also been assumed that both the optical and the atomic heat 
baths have zero thermal occupation at ooc and Ua, which is a good approximation at 
room temperature for Uc and uja within the visible-light frequency range. The population 
inversion of the atom is calculated from the density matrix via AP = Tr [pd^], with 
“Tr” being the trace of operators. 

We first check the validity of our analytic model for the sinusoidally modulated Rabi 
oscillation. In each of Figs. 2 a and 3a we plot three curves: Eqs. (IH]) [blue dash] and ([9]) 
[red dash-dot], and the numerical simulation of the master equation ffTOj) [black solid] 
with pit = 0) = |0)^ |e)^ |0)^^ (0| ^ (e| ^ (0| and k = 7 = p = 0. In Fig. 2a, we take the 
parameters to be coa = oje, gca = O.Swm, gem = O.ltUm, for which Eq. ([9]) is equivalent to 
the more general Eq. ([ 8 ]). We find that the two analytic curves well fit the numerical 
one, and the system undergoes a perfect sinusoidally-modulated Rabi oscillation. In 
Fig. 3a, we set Ua = oOc, gca = O.dScUm, Pem = O.lUm- For these parameters, Eq. ([9]) is 
no longer valid since ojm 7 ^ “^gca, but Eq. (jS]) still holds. In this case, Eq. (jS]) agrees 
well with the numerical result, and the modulation of the Rabi oscillation is not ideally 
sinusoidal. Fig. 2a and Fig. 3a conhrm that our analytic model correctly describes the 
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Figure 4. Modulated Rabi oscillation in presence of thermal effects. The parameters 
are gem = O.lUm, CUa = LOc - O.OluJm, gca = 0.49u;m, H = 0.02wm, 7 = 0.005wm, 
/i = 2 X nth = 10, and the initial thermal phonon number of the mechanical 

resonator is =0.5. 


coherent dynamics of the hybrid system within the proper parameter range. 

In Figs. 2b and 3b, we plot the numerical results for the temporal evolution of 


the phonon number = Tr 


pOb 


with the same parameters as in Figs. 2a and 3a, 
respectively. According to the analysis in the last paragraph of Sec. Ill, the Rabi- 
oscillation amplitude is maximum (minimum) when the phonon number is minimum 
(maximum). Comparing Fig. 2a and Fig. 2b, as well as Figs. 3a and 3b, immediately 
conhrms this conclusion, and thus provides further support for our analytic model in 
the previous section. 

Next we investigate the situation where dissipations are present by adopting 
K = 0.02um, 7 = O.OOSwm, /i = 2 X nth = 10. We also allow that the other 

parameters do not rigorously meet (but not deviate too much from) the requirements 
for Eq. ([9]), with oja = oJe — O.lcjm, Qca = O.dQcUm, and gem = O.lcUm- Moreover, since 
the mechanical mode can at best be pre-cooled to a small but hnite mean occupation 
number < nth, we set the mechanical resonator to be initially in a thermal state 
(rather than the ground state) of n[^ = 0.5. The numerical solution for AP, as plotted 
in Fig. 4, unambiguously demonstrates that the modulated Rabi oscillation survives 
even though it is siginhcantly degraded from the ideal one illustrated in Eq. ([9]) [or in 
Fig. 2a]. 


5. Summary 

In summary, the quantum behavior of an excited two-level atom in a vacuum 
optomechanical cavity can be qualitatively modihed by the mechanical resonator. 
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Particularly, if the atom-cavity Rabi splitting is on resonance with the mechanical 
mode, then the standard Rabi oscillation of the atom is sinusoidally modulated. The 
modulation originates from the “polariton-phonon” transition caused by the coupling 
between the optical held and the mechanical resonator. We explained this phenomenon 
with an analytic model in a three-dimensional Hilbert subspace of the hybrid system, 
which was further conhrmed by numerical simulations of the density-matrix master 
equation. The modulated Rabi oscillation was found to be reasonably tolerant on 
thermal effects and non-idealized parameters. 

6. Acknowledgement 

We gratefully acknowledge hnancial supports from the National Natural Science 
Foundation of China (No. 11375081, 11347026), the Shandong Provincial Natural 
Science Foundation (No. ZR2013AL007, ZR2013AM012), and the Start-up Fund of 
Liaocheng University. 

[1] S. M. Dutra, Cavity Quantum Electrodynamics: the Strange Theory of Light in a Box (John Wiley 

& Sons, Inc. 2005). 

[2] S. Haroche and J. M. Raimond, Exploring the Quantum: Atoms, Cavities, and Photons (Oxford 

University Press 2006). 

[3] I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, “Theory of ground state cooling of 

a mechanical oscillator using dynamical backaction,” Phys. Rev. Lett. 99, 093901 (2007). 

[4] F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, “Quantum Theory of Cavity-Assisted 

Sideband Cooling of Mechanical Motion,” Phys. Rev. Lett. 99, 093902 (2007) 

[5] T. J. Kippenberg and K. J. Vahala, “Cavity Optomechanics: Back-Action at the Mesoscale,” 

Science 321, 1172 (2008). 

[6] M. Aspelmeyer, P. Meystre, and K. Schwab, “Quantum optomechanics,” Phys. Today 65, 29 

( 2012 ). 

[7] M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, “Cavity optomechanics,” arXiv:1303.0733 

(2013). 

[8] F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity Optomechanics with a Bose-Einstein 

Condensate,” Science 322, 235 (2008). 

[9] P. Rabi, “Photon Blockade Effect in Optomechanical Systems,” Phys. Rev. Lett. 107, 063601 

( 2011 ). 

[10] A. Nunnenkamp, K. Borkje, and S. M. Girvin, “Single-Photon Optomechanics,” Phys. Rev. Lett. 

107, 063602 (2011). 

[11] K. Stannigel, P. Komar, S. J. M. Habraken, S. D. Bennett, M. D. Lukin, P. Zoller, and P. Rabi, 

“Optomechanical Quantum Information Processing with Photons and Phonons,” Phys. Rev. 
Lett. 109, 013603 (2012). 

[12] P. Komar, S. D. Bennett, K. Stannigel, S. J. M. Habraken, P. Rabi, P. Zoller, and M. D. Lukin, 

“Single-photon nonlinearities in two-mode optomechanics,” Phys. Rev. A 87, 013839 (2013). 

[13] T. T. Heikkila, F. Massel, J. Tuorila, R. Khan, and M. A. Sillanpaa, “Enhancing Optomechanical 

Coupling via the Josephson Effect,” Phys. Rev. Lett. 112, 203603 (2014). 

[14] Jie-Qiao Liao and Franco Nori, “Single-photon quadratic optomechanics,” Scientific Reports 4, 

6302 (2014). 

[15] K. Hammerer, M. Wallquist, C. Genes, M. Ludwig, F. Marquardt, P. Treutlein, P. Zoller, J. Ye, 

and H. J. Kimble, “Strong Coupling of a Mechanical Oscillator and a Single Atom,” Phys. Rev. 
Lett. 103, 063005 (2009). 


Sinusoidally Modulated Vacuum Rabi Oscillation... 


10 


[16] J. Restrepo, C. Ciuti, and 1. Favero, “Single-Polariton Optomechanics,” Phys. Rev. Lett. 112, 

013601 (2014). 

[17] O. Kyriienko, T. C. H. Liew, and 1. A. Shelykh, “Optomechanics with Cavity Polaritons: 

Dissipative Coupling and Unconventional Bistability,” Phys. Rev. Lett. 112, 076402 (2014). 

[18] J. M. Pirkkalainen, S. U. Cho, F. Massel, J. Tuorila, T. T. Heikkila, P. J. Hakonen, and 

M. A. Sillanpaa, “Cavity optomechanics mediated by a quantum two-level system,” Nature 
Communications 6, 6981 (2015) 

[19] M. Scully and M. Zubairy, Quantum Optic.s (Cambridge University Press 1997). 



